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A new set of projection operators is constructed to suitably handle non-relativistic theories of 
^SJ , gravity with anisotropic scaling, including the ones with parity-violating terms. This alternative 

procedure allows us to discuss unitarity and spectral properties for different formulations of the 

■ Hoi'ava-Lifshitz gravity. This task, that generally involves lengthy algebraic steps, becomes more 
04 \ systematic and greatly simplified in terms of the projectors we work out. Moreover, this procedure 

5^ . allows us to fix the number of propagating degrees of freedom and the structure of gauge symmetries 

■ is readily determined. In order to test the efficacy of the technique at hand, the unitarity and low- 
energy regime of a general Hofava-Lifshitz Gravity model are investigated. 
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I. INTRODUCTION 

(~| \ The recent proposal of Petr Hofava for a new theory of gravity [ll-0| was received with great enthusiasm by the 
Field Theory, Gravitation and Cosmology communities [5l-l39f and interesting consequences for cosmological models 



^ ] have been worked out |40l - |47| . The large appraisal for this model, dubbed as Hofava-Lifshitz gravity (HLG), came 
^ . from the fact that it is a candidate for a consistent quantum field theory for gravity, since it has a remarkably 
(y-j ' improved ultraviolet behavior. However, to achieve such a result, Hofava had to give up diffeomorphism invariance 
, as a fundamental symmetry so to impose anisotropic scaling of the space and time dimensions. 

■ Despite its success, it soon became evident that Hof ava's orig inal formulation including the detailed balance condition 
l/^ [ is not obviously a phenomenologically viable model (Tsl. [Sll. [32| . So, to reconcile HLG with experimentation, it became 
• necessary to loose up such condition and allow the inclusion of a large number of Lagrangian terms with arbitrary 
coefficients in the model or introduce auxiliary fields to enlarge the symmetries of the model. Besides this, HLG can 
be formulated in two distinct versions: nonprojectable or projectable - which refers to the dependence or not of the 
lapse function with the space coordinates. These different possible formulations brought a discussion in the literature 
concerning whether each of these formulations should be considered as a good, bad, healthy, or ugly version of HLG 

rS , All this debate included not only analysis of the phenomenological aspect, but also considerations on the the 
^ ■ quantum consistency of each version of theory. It is in this context that we judge necessary (and convenient) to 

" " ' develop a general method for the attainment of the propagators and, thereby, the identification of the excitations 
present in the possible models based on the HLG frame. With the propagators at hand it is possible to systematically 
obtain a description of the particle spectrum and the relativistic and quantum properties of scattering processes. More 
exactly, one is able to pinch out the propagating modes and determine conditions on the coefficients of the Lagrangian 
so as to impose the absence of tachyons and ghosts. 

There are several methods for the attainment of propagators, but, particularly in the case of weak field approxima- 
tion for quantum gravity, which is our interes t, algeb raic methods have been intensively developed, specially the one 
based on the spin projection operators (SPO) |48l - [5(]j . The SPO for Lorentz invariant models has the interesting prop- 
erty of being a orthonormal basis that decomposes the fields into definite spin-parity sectors. In this paper we follow 
the lines of previous works [sil - ls^ l to build up an orthonormal basis of operators that is specially suitable to handle of 
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HLG models, including the ones with parity-breaking terms. With this result we may discuss the spectral properties 
of the HLG and understand better the extra 3D spin-0 propagating mode that appear in the theory. Moreover, this 
technology allows us to determine the number of propagating degrees of freedom and the gauge symmetry structure is 
easily determined. We would like to point out that recently F. S. Bemfica and M. Gomes also discuss the propagators 
and consistency conditions on the spectrum for special classes of HLG models by adopting an independent method 

This paper is organized as follows: In Sec. |TT1 we build up the basis of projection operators suitable with the 
symmetries of HLG. We present in the Sec. IHIl a general method for the consistency analysis of the particle spectrum. 
The Lagrangian model of our interest shall be described in the Sec. IIVI and the attainment of its propagators along 
with the spectral analysis is done is Sec. |Vl We present our Concluding Assessment in Sec. IVII emphasizing in some 
consequences of the propagator structure of HLG in the low energy regime. In the Appendix [Xj we list the inverse 
of the spin-0 matrix used to calculate the propagators of Sec. |Vl whereas the Appendix [B] is aimed to derive some 
useful algebraic relations that the projection operators satisfy. 

Unless stated otherwise, we shall use Latin letters (a, b,c, . . .) for 3D spatial index. We also define Sab = diag (1, 1, 1) 
and eabc as the completely antisymmetric symbol normalized to eoi2 = +1- Also, we shall follow the conventions: 
R'bcd = dcTtd + ^ce^d - (c ^ rf), = ig'"* {dbgdc + d^gdb - dagbc), Rbd = R%ad^ R = g'^'Rbd, with gab being the 
spatial metric. 



II. DECOMPOSING THE GRAVITON EXCITATION MODES 



In the Hofava-Lifshitz gravity (HLG), the breaking of (1 + 3)-dimensional diffeomorphisms invariance is imple- 
mented by endowing space-time with a preferred foliation of three-dimensional spacelike surfaces. This defines the 
splitting of coordinates into space and time that explicitly breaks general covariance down to the subgroup of coordi- 
nate transformations 

t^i{t), x'' ^ {t,x''), (1) 

referred in the literature as foliation-preserving diffeomorphisms. 

To fulfill the task of defining a Lagrangian model with such a symmetry, it is natural to make use of the ADM 
(Arnowitt-Deser-Misner) decomposition, which in the convention we are adopting, decomposes the space-time metric 
'"^^ gpif ill terms of the lapse A/", the shift Na, and the spatial metric gab as 

ds^ = (TV^ - NaN") dt^ - 2Nadx''dt - gabdx^dx^. (2) 

From the equation ([2]) we may relate the inverse of the four-dimensional metric and the ADM fields by 

(4) 00 _ J_ (4) 0a__^ W ab _ ab , (o-. 



and V^detW^ = Afy/det gab- 

Having defined the fundamental fields of the model, we may perform a weak field expansion around a Minkowski 
metric solution in order to study the graviton excitation modes of HLG, 

TV = 1 + n, gab^ Sab + Kb, (4) 

with n and hab being the quantum fiuctuation of the lapse function and the spatial metric, respectively, and the shift 
function Na being the quantum fluctuation itself. 

Breaking down the Minkowski symmetry, 5*0 (1, 3), of the linearized theory to a SO (3) has some nontrivial con- 
sequences, specially in the interpretation of the propagating modes. An usual relativistic quantum field theory is 
invariant under Poincare transformations. It provides a classification for particles as unitary irreducible represen- 
tations of the Poincare group by the quantum numbers: mass and spin. In particular, the spin is characterized by 
the unitary representations of the 5'0 (1,3) little group, that is the subgroup of the Lorentz group that leaves the 
representative four-momentum unchanged. 

For the class of HLG models, we may proceed in a strict analogy. We split our propagating modes in "what would be 
the spin" defined the subgroup of the rotation group which leaves the representative three-momentum invariant. So, 
to arrange the degree-of-freedom into these "spin modes", we may develop a complete set of projection operators for 
Hoi'ava-like models. For this task, it is reasonable to make use of the already developed three-dimensional SO (1,2) 
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Barnes-Rivers projection operators \5)\ and adapt them for a SO (3)-invariant model implementing the following 
substitution: 



VfiC' ^ Sab J y ka 



(5) 



where /z, j> = 0, 1, 2, a,b = 1, 2, 3, fc/i is the (l+2)-D relativistic momentum and ka is the 3D spatial momentum. For 
example, in the given situation, the transverse {9ab) and the longitudinal {cjab) operators, which are building blocks 
operators become 



Pah 
<jJab 



— Sab — OJab, 
kakb 



(6a) 
(6b) 



with fc^ = kaka- 

In this way, we may organize the HLG projection operators in a matrix-form according to the respective spin sector 
(P (0), P (1), and P (2)) to cover the case of parity- preserving Lagrangians containing the fields n, Na, and hab- They 
are arranged by the possible spin representations carried in the quantum fields - n G 0, Na € ® 1, hab € © 1 © 2, 
which comes from usual group-theoretic argumentation. So, we cast them as 



P(0) 



hab 
hab 
Na 

n 



had had Na _ n 

2^ab0cd ■:y^(^ab^cd -^^abkc ab 

-^^obQcd 

'y^kaObc 



^ab^cd 


^abkc 


^ab 


kai^bc 


^ab 


ka 


^cd 


h 


1 



(7a) 



P(l) = 



1 
%/2 



had 

bc^ad 

kn 



Na 



{Qac^bd + dbc^ad + dad^bc - 
- Ocah) 



Sbd^ac) {Oackb + ^'bc^a 



Tab 



(7b) 



^ " (2) — 7: [Oac^bd + dadObc — Oab^cd) 



(7c) 



where we have defined ka = 

However, even in Hofava's original proposal with detailed balance and its subsequent modifications there appear 
parity violating terms, that notably cannot be written in terms of the SO (3) Barnes-Rivers operators. A hint to 
circumvent this difficulty is to note that in the presence of parity violating terms, the transverse modes may propagate 
independently and not as parity-doublets (as discussed in Ref. [s^l and put into test in Ref. (Isj). In practice, the 
two degrees of freedom of the ^-operator can be split by two orthonormal operators. 



9a& — Pab + CTab 



(8) 



satisfying = p, — a. and p • tr = 0. This decomposition is useful to handle with parity violating term such as 



the ones containing the e-symbol, since p and <j may be defined to satisfy 

e'^^^fcc = i^{pab - (Jab) ■ 



(9) 



The degree-of-freedom projection operators for relativistic planar models {SO (1, 2)-invariant theories) with parity- 
breaking terms have been thoroughly worked out in the Ref. [52] and can be used adequately for HLG models with 
parity violating terms. Specializing these projection operators for the fields of the ADM decomposition, one obtains. 



P(l) 



hab ( + ) 
hab (-) 
Nai + ) 
Na{~) 



had{ + ) had{-) _ Na( + ) 

2pacUJbd _ 2eghePa'Jci^bdk'' \/2packb 

2egheCraPc'^dbk'' 2aaci^bd V2i 

V2pbakc V2eghe<7bPli^ec Pab 



V2e 



Na{- 



V2( 



^ghePbgO'a'^ce 



ieghe<j'f!;ifa^ec 

V2abakc 



CgheP 

f^gheCTaPc^eb \/2a aj^b _ 
(.^B^PafOhgkh 



-e'^^^CJadPbekf 



Cab 



(10a) 
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P(2) = ('"''S+J 



hcd(+) h^ai-) 
^{PadPbc + (^adC^bc - Pab<^cd — labPcd) <^ghe{PcaCrd Pb ~ cb Pd'^a) 
eghe{PacO-b P^a ^ <^adPl(7c W '^PacCTbd 



(10b) 



In the above matrices it is understood the operators share the same symmetrization properties (with numerical factor) 
of the associated fields. 

It is noteworthy to emphasize that all the HLG projection operators ([7a|l - (|7c|) and (|10ap - (jlOb[) have been constructed 
in such a way that the orthogonality relationship and decomposition of unity relations are fulfilled: 



^ P^J{I)c.f^Pkl{J)f}-, = S,kPu{I)c.^, /, J = 0, 1, 2 

Pii{J)alS — Sal3, 



(11a) 
(lib) 



i,J 



where the Greek indices a, /3, 7, . . . represents here either tensorial, vector or no spatial indices. 



III. ANALYSING THE CONSISTENCY OF THE PARTICLE SPECTRUM FOR A GENERAL MODEL 



To understand the convenience of orthonormal basis of operators for the attainment of the propagators, let us work 
out a general spectral analysis for models suitable with the symmetries defined by HLG models. For such, let £2 be 
a quadratic Lagrangian for a HLG model: 



C2 ^^iPaOa/Sipp +^'PaSa, (12) 
Q/3 a 

where the wave operator Oap is a local differential operator, tpa = {n, Na, hab} carry the fundamental quantum 
fields of the model, and Sa = {S,Sa,Sab} are the matter sources of the lapse, shift, and spatial metric, respectively. 
We can systematically analyze the spectrum and unitarity of this model by means of a decomposition of the wave 
operator Oap in terms of the HLG projection operators P^^ (J) in the momentum space: 



£2= E (A, (13) 

where a {J)ij are the coefficients of the expansion of the wave operator. The diagonal operators P^"^ (J) are projectors 

in the spin- J sector of the field (p and the off-diagonal operators P^^ (J) (i 7^ j) implement mappings inside the spin 
subspace among the fields (p and ^p. 

The properties (|llap - (jllb[) of the HLG projection operators reduces the task of inverting the wave operator into 
a straightforward algebraic exercise. Indeed, if the coefficient matrices a{J)ij are invertible, then the propagator 
saturated with the external sources Sa is given by 



n = z ^ 5>(J),/P„-(J)„^5^. (14) 

Nevertheless, the coefficient matrices a iJ)ij may become degenerate due to a gauge symmetry of the Lagrangian 
model. This poses no serious problem to the attainment of the propagator since, in a such case, there are sources 
constraints that consistently appear in order to inhibit the propagation of non-physical modes. Both, the gauge 
transformations of the fields and the source constraints may be obtained from the degeneracy structures of the 
coefficient matrices, 

Scba = ^^^""'"^ (J) , for any j (15a) 

J2 ^/^'"^ (J) P^3 iJ)af! S'p (J) - 0, for every i and J (15b) 

/3j 
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with fp (J) being arbitrary functions and being the right (left) null eigenvectors of a {J)ij- 

With it in mind, one should be convinced that the correct procedure for the attainment of the propagator is to invert 
any largest nondegenerate submatrix of a {J)ij- In practice, it suffices to delete the degenerate rows and columns of 
a{J)^j (which we denote A(J)^^) according to the number of gauge symmetries, to invert A(J)j^, and saturate it 
with physical sources, in order to obtain the propagator: 



n = z ^ s:AiJ);/p,,iJ)^^s'fs. (16) 

The propagators provides key ingredients for the understanding the scattering processes of the theory under in- 
vestigation. In particular, we may analyze more carefully some conditions to guarantee the tree-level unitarity and 
causality of the Lagrangian model. In general, the poles of the propagators relation for HLG model is 

w2_Q(fc) = 0, (17) 

where Q may be a polynomial or rational function of the absolute value of the spatial momentum k. So, for each 
propagating pole we will require that Q (fc) > 0, for every momentum k, to ensure it is not a tachyon mode. However, 
at low energies the higher momentum terms becomes less relevant and the dispersion relation takes a simple form 

^2 - c2fc2 = 0. (18) 

So, for each propagating pole we will require that > to ensure it is not a tachyon mode. 

On the other hand, the condition imposed by quantum field theory for the absence of propagating ghosts reads 



9mRes(n|po,e) > 0. (19) 

Fortunately, with the projection operators formalism we can take advantage from the general decomposition of the 
degree-of-freedom HLG projection operators, 

P^jiJ)c.p^(.-^r^^a^i^f, (20) 

where p is the parity related to the spin operator, to rewrite the propagator (|14p as 

H = ^i~^fJ2^^^(J^Q)^/^'A^'-Qi'')y'^ (21) 

where Sj = ip'£'S°' and A (J, Q)"> are the inverse sub-matrices with the pole extracted. Therefore, the positiveness 

condition for arbitrary sources, is ensured by the positiveness of the eigenvalues of the A{J^Q)~^^ matrix. 

Nevertheless, it can be shown, for each case, that these matrices has only one non-vanishing eigenvalue at the pole, 
which is equal to the trace of A {J, Q)^^ Therefore, the condition for absence of ghosts for each spin is reduced to: 



{-\YiYA{J,Q)-l\poie>Q- (22) 

This technical result converts the task of probing the unitarity, that is in general a time-consuming work, into a 
straightforward algebraic exercise of analyzing the matrices of coefficients. 



IV. DESCRIPTION OF THE LAGRANGIAN MODEL 



Having discussed a general method to obtain the propagator and to analyze the spectral consistency, we may focus 
our efforts in specific HLG models of our interest. Since one can impose additional symmetries and restrictions to the 
theory it must be made clear which version of the theory one is referring to. In Hoi^ava's original proposal, probably 
concerned in maintaining its simplicity, he mainly discussed the version with the detailed balance condition. This 
condition imposes that the possible spatial metric-dependent terms are of the form Sv — J dtd^ xAf ^^"'''QabcdS'''^ , 
with itself following from a variational principle ^S""^ — ^^^V^- Therefore, it implies in a huge reduction in 
the arbitrary parameters of the model, however it has been shown it becomes very difficult to reconcile this condition 
with phenomenological considerations 0, [H, HH H^l ■ 
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Other requirement that significantly alters the structure of HLG models is the projectability condition. It enforces 
the lapse field to be just a function of time Af = {t). Evidently, this avoids Lagrangian terms which are spatial 
derivatives of J\f, and thus simplifying the theory. A motivation for imposing this assumption, comes from the 
possibility of setting A/" = 1 by a gauge fixing. In the following, we shall compare how this condition affects the 
spectral properties of the model, but without the introduction of auxiliary fields. 

Although interesting considerations can be made with both conditions, we shall be mostly concerned with the 
version where the detailed balance and projectable conditions are neglected. Allowing numerous possible Lagrangian 
terms may turn the attainment of the propagator a painful task, but with the HLG projection operators this can be 
done with a systematic and simplified technique. This said, let us enumerate the possible Lagrangian terms compatible 
with the foliation-preserving diffeomorphisms (Eq. (U}). 

First, we start with KabK"-^ and K'^ = {^g°-''Kabf , which are written in terms of the extrinsic curvature 

Kab = ^ [gab " VaiVfc - VfeiVa) . (23) 

Kab is covariant under spatial diffeomorphisms and transforms as a scalar under time reparametrization. Here and 
throughout the manuscript, the dot refers to time derivative and Va is the covariant derivative associated with the 
three-dimensional metric gab- KabK""^ and = (^g°'^Kab) are referred kinetic terms, since they contain (at most 
two) time derivatives. 

The potential terms do not contain time derivatives. The ones built with the spatial metric gab have been enumerated 
in Ref. (3l| . Since in three-dimensions, the Riemann tensor is equivalent to the Ricci tensor the independent possible 
terms up to dimension [nf are i?^, RR^R^, R\R\R''a, RV^R, V^R, Vai^fccV"i?''^ e'''"=RadDbR/, i?^ R''''Rab, 
V^i?, e-^^- (rtdbT^^, + IrLrir^.V R, and 1. 



cd 3^ ad^ be^ cf 

However, the terms above do not exhaust all the possibilities. In fact, an important novelty brought by Bias, Pujolas, 
and Sibiryakov in Refs. d, Q was realizing that terms involving the object Aa = M~^daJ^ should be included since 
it is invariant under foliation-preserving diffeomorphism. The terms mixing Aa with the extrinsic curvature Kab with 
dimension [n] or less are K°-^AaAb, K^^S/aAb, and KVaA". They have the property to be odd under T and CPT 
transformations. And the list of potential terms built with Aa reads AaA"^, Aa^'^A'', {AaA'^f, AaAbR''^ Aa^'^A" , 
{AaA"')'^ , AaA'^AbAcR'"^--- Notably, with this improvement it has been argued that HLG does not run into problems 
with strong coupling and could be a phenomenologically viable model [1, 0|. 

After this brief consideration on the development of HLGs (for a deeper discussion we address the reader to Refs. 
0, 0-i, [11 [13, m, nil), we set up a general HLG model only with the terms that affect the propagator, 

Shl = J dtd^x^N {aKabK"'' + PK^ + Cy + Ca) , (24) 
where the potential terms related only to the spatial metric are 



Cg = aiR + biRabR"^ + ciR^ + ail^R + b2RabAR''^ + C2RAR + agA^i? + fiCcs + XCrc, (25a) 
where the Chern-Simons like term reads 

^cs = e^'^ (rtdbTi, + ^r{,r^,r^^.) (25b) 

and the Ricci- Cotton term is cast as 

Crc = e'^^'^RadDbRf. (25c) 
The Lagrangian with the terms related with Aa reads, 

Ca = rK''^VaAb + sKVaA"- + uiRVaA°' + u2ARVaA'' (25d) 
+ riiAaA" + ri2AaAA"' + 

We should emphasize that, at this point, we shall adopt the posture of not avoiding the discussion of the terms 
that violate P, T or CPT symmetry. Up until now, we have no fundamental reason for leaving them out and it is our 
interest to understand how they affect the propagating modes dispersion relations. 
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V. ATTAINMENT OF THE PROPAGATORS AND SPECTRAL ANALYSIS FOR THE GENERAL 

HORAVA-LIFSHITZ GRAVITY MODEL 



Having established the model of our interest, we may work in the wave operator defined by Eq. ([24|) . After some 
algebraic calculations (including partial integrations), it is possible to bring the quadratic Lagrangian for the action 
(|24| into the form of the Eq. ((T2|) . With this procedure we determine the corresponding wave operator. In the sequel, 
we perform the wave operator decomposition in terms of the HLG projection operators, as indicated by Eq. ()13|) . 
aided by the relations developed in the Appendix |B] 

Some remarks are in order. In the attainment of the spin-1 matrix of coefficients one realizes that there is no need 
for decomposing the propagating modes its degree-of-freedom, since it appears no parity violating contribution in 
this sector. For simplicity, we write them in terms of the parity- preserving operators (j7bp . On the other hand, the 
Chern-Simons and Ricci-Cotton parity breaking terms (Eqs. (|25bp and (|25cp . respectively) imply the violation of 
parity in the spin-2 sector. It forces us to use the degree-of-freedom operators (llObp in this sector. With this in mind, 
we write out the coefficients matrices of the wave operator expansion. 



a (2) 



2 



he) e 



(26a) 



a(l) 



-aijj- 



V2 



V2 



ae 



(26b) 



/ (f + 



a(0) 



13) w2 + lAk^ 



1^ I ,2 



i (a + /?) 



■~l3ojV2k^ V2 {a - uP) k^ + ^siujk^ \ 
{a + fi) LuVe i (r + s) iuk'^ 



-l3ujV2k^ -{a + 13) LuVe 2 {a + 13) ^ 

\ V2 {a - uP) e - ^siLje -^{r + s)iujk^ {r + s)ik'^Vk^ 



(r + s) ik^ve 



(26c) 



where we have redefined some coefficients to avoid cluttering the notation: 



(27a) 
(27b) 
(27c) 
(27d) 
(27e) 
(27f) 

As mentioned in Sec. |lTl the degeneracy structure of the matrices (|26ap - (|26cp brings us information about the gauge 
symmetries and the source constraints of the model. In particular, the spin-1 matrix is degenerate and the 2"** and 
3'"'* rows of the spin-0 matrix are proportional. Adding together the result of equation (jl5b[) for these two matrices, 
one obtains 



a = 


ai 


- a2fc^ 




b = 


bi 






c — 


Cl 


- C2fc^, 




u = 


Ul 


— U2k^, 




V = 


Vi 






A = 


a-\ 


- (36 + f 


Ic) e. 



kaSab = 0. 



(28) 



and all 



For a consistency test, one may verify that for the Einstein-Hilbert gravity (EHG) (a = — /3 = ai — i^^q 
other coefficients vanishing) there appears another symmetry in the spin-0 matrix: the 3'"'' and 4*'' columns and rows 
will become proportionate. It will give rise to another constraint among the sources 



ujS + kaSa = (only for E-H). 



(29) 



These results are compatible with the relativistic source constraint k^S^'^ = (/x, = 0, 1, 2, 3 ) if one takes into 
account that: (i) the contributions of /i"" and are encoded in N"': it generates a factor of 2 in the matrix 
decomposition, which implies in a ^-factor in the source iSq and (ii) the weak field expansion M — 1 + g, along with 
the definition in the ADM decomposition of implies that is identified with 2n; this also brings a i-factor 
associated with the source S. 

With these matrices coefficients, we may follow the analysis of the spectral consistency as proposed in section |lTl 
For the sake of clarity, we separated the discussion for each spin mode. 
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A. Spin-2 Sector 



Let us analyze is the spin-2 mode of the graviton. The matrix a (2) from equation p6ap is non-degenerate. To 
obtain the propagator (|14p its inverse may be readily obtained 



a 



where the denominator D(2) = (f'^^ ~ 5 ("^ ~ ^^^) ^'^Y ~ (m+ h'^k^Y ^^^^ ^wo poles in the energy. This is an 
interesting feature of the parity breaking terms: they allow each spin component to propagate independently, with 
different dispersion relations 

ujI = - (a-bk'^)k'^ ±- ( fi + -XkA k'^Vk^. (31) 
a ' a \ 2 y 

Using the posture of requiring a positive speed for every propagating mode, one must ensure that 

C(2) = - > 0. (32) 

The condition for the absence of ghosts must be made for each pole utj. = Q± (k). In either cases, implies that 
1/a > 0, so the ghost and tachyon free conditions are given by 

Spin-2 : oi > 0; a > 0. (33) 

From this analysis we conclude that the Chern-Simons and Ricci-Cotton terms, which are related with the pa- 
rameters fi and A, do not interfere in the unitary relations as its related terms in the three-dimensional topological 
gravity. 



B. Spin-1 Sector 

Since the spin-1 matrix is degenerate, we shall use the prescription of the equation (jl6p . A largest non-degenerate 
matrix can be obtained by deleting its 2"'' columns and rows. Its inverse is given by 



This mode poses us with a strange pole o;^ = 0. However, is related to the projection operator Pl''^^ (1) = 
i {Oac^^bd + (^bc^ad + (^ad^bc + dbd^ac): which implies that the spatial momentum ka is contracted with the source 5°^ 
in the expression ([TC)) . In this way, using the source constraint one verifies that the residue at the pole of 

the propagator saturated with physical sources vanishes identically. We must therefore understand this pole as a 
non-propagating mode. 



C. Spin-0 Sector 

We left for last the spin-0 sector, which brings the most interesting discussions in the literature. It is a natural 
consequence of the lost of symmetry that new propagating modes appear and this is reflected in this sector. However, 
it is of a major importance to verify the spectral consistency of this mode. In fact, the EHG augmented with 
relativistic higher-derivative corrections has been shown to be renormalizable fs^, but exactly the 4D spin mode 
yields a ghost excitation which jeopardizes unitarity of S-matrix and its possibility to become a consistent quantum 
model for gravitation. 

The second and third columns of the matrix a (0) (equation (|26cp ) are degenerate, due to a residual gauge symmetry. 
The inverse of A (0) is cast in the Appendix |A] and stems a dispersion relation. 



1 


(iAj7- (a-uP)') 


1 (a + /?) - \AP (r + sf 


ia?7 (a + 3/3) - i (a ((r -f sf + 2s2^ 


1 + 2/3r2) 


1 fc2 



0. 



(35) 
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This great number of arbitrary parameters makes the spectral analysis a complicated tasks and many intriguing 
situation occurs. In the special case where 7] = u = 0, and r = — s 7^ we obtain the dispersion relation u'^ + 4^- = 0, 
which is a constant negative energy solution. 

If the CPT violating terms are absent r = s = 0, we have no great difficulty to proceed. In this situation, 



/ 2r]{a + /3) -2^2?]^ -V2a (a + /3) 

A(Or = -^ -2v/2ry/3 M^+m^^-{^a^-2A,y ^a/S |, (36) 

with D(o) — ri{a + 3/3) aui^ — (2a^ — Ar\) [a. + ji) and the dispersion relation is given by 



2 {2a^-A7^) (a + /3) 



- (\or{ k'=0, (37) 

7] {a + Sp)a 

which in general is a rational function of the spatial momentum fc^ . It recovers the dispersion relation obtained in the 
reference Q. 

For small values of momentum the positivity of the propagation speed for this mode reads, 



m {a + 3/3) a ^ ' 



To obtain the ghost-free condition, one ought to verify the residue of the matrix at the pole. Replacing the dispersion 
relation (j37p and after a inspection (in a CAS, for example) one verifies that A (0) ^ \poie has only one non-vanishing 
eigenvalue, which is equal to its trace, so the expression (I22p is valid, 

t 4,0^-1, (a^ + 2.y^)(a + /3)V4/3V 

tr A O) po/e = 2/ ,R\f I QflN • (39 

arj^ [a + p) [a + Sp) 

111 this way, the condition for absence of ghosts can be stated as 

a{a + P){a + 3^) > 0. (40) 

Combining the results ((38|) and (|40p with the spin-2 condition ([33l) one must impose that 

Spin : {a + (3) {a + 3/3) > 0; 771 (2ai - rji) > 0, (41) 

for the absence of ghost and tachyon in this sector. 

For comparison, let us come back and analyze how the projectability condition affects the matrix structures and 
thus the propagating modes. Imposing the restriction that the lapse field does not depend on the space variables 
J\f = JV (t) implies that the terms in the quadratic Lagrangian containing the perturbation of the lapse n are total 
derivatives. Thus, the fourth column and row of the spin coefficient matrix vanishes identically, 



a{0) 



i(a + /?)w2 -{a + /3)ioV¥ 



V 



-/3a;V2fc2 




{a + l3)ujVk^ 2(a + /3)fc2 
/ 



(42) 



With this simplification, one readily realizes that a (0)^ is doubly degenerate: the column 4 is null and the columns 
2 and 3 are proportional. A largest non-degenerate matrix of a (0)^ is obtained by deleting the 3'"'' and 4*'' columns 
and rows, 



^(0), 



>^ 1 ,2 



i (a + /?) w2 



Its inverse is given by 



1^ , ,2 
2 ' ^AP 



(43) 



(44) 
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with 



D(o)p = (a + 3/3) auj'^ + (a + /3) Ak^ 



(45) 



providing the dispersion relation for this propagating mode. To ensure that this mode is not a tachyonic excitation 
we must impose 



-(0)p 



(a + 13) ai 



> 0. 



(46) 



It can be explicitly verified that Res |a (0) ^| 



pole 



{a + 3/3) a 

is degenerate, in such a way that its only non- vanishing eigen- 



value matches its trace. So one is enforced to impose that 



tr (^Res|A(0)"^} 

The only way to reconcile the conditions and (|47p would be to impose ai < 0; but, this contradicts the unitarity 
condition for the spin 2 sector. Therefore, we may infer that the HLG with the projectable condition is not compatible 
with perturbative unitarity at the tree-level approximation. 

An alternative to circumvent this issue is to enhance the model with an Abelian gauge symmetry 0, [T^ [soj , in 
order to inhibit the propagation of the spin-O mode . This is possible with the introduction of a gauge field A and an 
auxiliary scalar field i/, known as the Newtonian prepotential. This approach yields a model with the same number 
of degrees-of-freedom as the EHG. 



nol.p. I n. ( n -\- B) ( cy 4- 38) 



> 0. 



(47) 



VI. CONCLUDING ASSESSMENT 



In this manuscript, we have built a set of orthonormal projection operators suitable for non-relativistic theories of 
gravity, including the ones with parity violating terms. The advantage of this construction for the attainment of the 
propagator and the spectral analysis, which is in general a time-consuming task, is that it can be done in a systematic 
and simplified way, including the handling of the gauge symmetries of the model. With this methodology, we were 
able to determine conditions over the coefficients of general Hofava-Lifshitz gravity (HLG) in order to restrict the 
propagation of tachyons and ghosts. 

With the general result of the propagators (Eqs. (1501) . ([M)) . and we may contemplate some results in a 

low-energy regime of the theory, where the terms depending on higher spatial momentum become less relevant. In 
such situation, the general action of the nonprojectable HLG (Eq. (j24p l may be simplified to 

Shlg = J dtd^x^N {aKabK"''' + f3K'^ + aR + TjAaA"') . (48) 
and the inverse of the coefficient matrices of the propagator for this simplified model are cast as 



2r]{a + f3) -2V2rjl3 -V2a{a + l3) 

AiOrHLG = J^\ 2r,i.+2,W-M2a-,)k^ 2a/3 I, (49a) 



^(D.L-i^, «(2)-,.-^, (49b) 

with L»(o) ='q{a + 3/3) aw^ - a {2a- 77) {a + I3)P. 

It is instructive to compare the propagators of the HLG (Eqs. (I49a[) and (|49bp ) with the Einstein-Hilbert gravity 
(EHG) propagators in ADM coordinates with HLG projection operators. They may be written as 
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-1 2 _i 



By a quick glance at these matrices of coefficients one may set up comparisons between both models. Analyzing 
the pole structure of the propagators, we may conclude that the spin-2 massless mode propagate in both models, but 
with different speed of propagation. This corresponds to the usual helicity-2 graviton mode of the relativistic model. 
Both models coincide by not having dynamical poles in the spin-1 sector. However, HLG has a propagating pole the 
spin-0 sector, which is absent in the EHG. 

This new degree of freedom may lead to odd effects. For example, the helicity-0 mode of massless limit of Pauli-Fierz 
gravity (PFG) generates a factor | in the bending of the light by the Sun, if one scales the coupling constant of both 
PFG and EHG such that the Newtonian limit holds |56l460( (see Ref. [61| for a recent review on the subject). This 
effect can be explained by the different way that the helicity-2 and helicity-0 modes couple to scalar and vectorial 
matter sources. One should be aware that this apparent paradox, known as vDVZ discontinuity, appears only in the 
perturbative approach of the linearized theory. 

Now, we are ready to compare the HLG and EHG in a situation of a point mass static source coupled to relativistic 
matter sources. The static point-like source can be described by 5 = MS{uj), Sa = Sab — 0. The propagator 
amplitude H for this source coupled to arbitrary sources {S' ,S'^,S'^f^} is 

- ^UHLG = MS iu) j^^^, { {Oab - J^fa?j 5^, - ^^j , (51a) 

-illEHG = MS (c.) ^ f^OabSU - , (51b) 

respectively for HLG and EHG. In the particular case of the Newtonian test with a point- like source of mass M' 
separated by the distance of i?, the energy of interaction is given by 

_ 1 MM' _ 1 MM' 

^^^^-"2(2a-r,) R ' ^^^^ " ^4^^' ^^^^ 

respectively for HLG and EHG. 

On the other hand, the relativistic electromagnetic matter should be dealt more carefully than just a scalar matter, 
for which 5^ = 0. In terms of the ADM-components the traceleness of the electromagnetic tensor reads, 

2^EM + ^abS^ti = 0. (53) 

Using the identity (1551) together with the source conservation constraint one can show that the interaction energy 
of the point mass with the light beam exceeds the value of Eq. (|52| by a factor of two for both theories (replacing 
mass with total energy). This implies that these experiments with static sources cannot distinguish the EHG from 
the nonprojectable HLG at the tree-level order of perturbation, since both are scaled by multiplicative constant. 

In this vein, one may look foward for a dynamical experiment that may unveil the effects of the extra scalar mode. 
An interesting possibility is to set up the test of gravitational pulsar experiment. It consists of a spherically symmetric 
harmonically pulsating mass distribution. In fact, according to the Birkhoff theorem, the EHG radiation from any 
spherically symmetric source is zero [6^]. Looking to the propagators coefficients of the EHG (Eqs. (|50ap - (l50bp ) one 
may verify that it is indeed the case, since the propagating pole of the EHG is associated with a transverse operator 
P (2) that vanishes when contracted with a spherically symmetric source. The propagator structure of HLG leads 
to a different situation, because it has non-vanishing coefficients with propagating pole associated with longitudinal 
operators. In this case only the spin-0 mode carries radiating energy and, hence, the Birkhoff theorem is not valid in 
the nonprojectable HLG. 

Our efforts here do not contemplate the incorporation of the cosmological constant in connection with the spectral 
analysis of the modes present in HLG. We understand that this is a fairly interesting issue, which may motivate further 
investigation. Indeed, in the recent paper of Ref. [26], the authors pursue a stimulating discussion on the stability 
of HLG black holes in connection with the AdS/CFT correspondence and possible applications of the holographic 
principle. In view of that, we shall be going further to endeavor a study of the consequences of the cosmological 
constant in the frame of the approach we report here. 

Also, the coupling of Lifshitz fermions to HLG could be an interesting matter in connection with our present 
investigation, since torsion degrees of freedom may be excited and their incorporation shall yield new interesting 
technical aspects in our treatment. The extension of our method to include torsion in the HLG frame is also motivated 
by the need to study the possible production and the consequent decay of massive gravitons at the TeV scale to finally 
compare with the results worked out in previous Lorentz-preserving quantum gravity models (63l-[69j . 
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Appendix A: Inverse of A (0) Matrix 

The inverse of A (0) matrix that is obtained by deleting the 3'''' row and column of a (0) (equation (|26cp ) is given 
by: 




(Al) 



where 



Dn 



A 



(0) 

11 



A' 



12 
(0) 



^13 



l(0) 



A 



23 
(0) 



^33 



^arj {a + 3/3) - ^k^ (a ((r + sf + 2s^^ + 2/3r^^ 



-Ar] ~ (a- uk^y ]{a + l3)- -Ak^ (r + sf 
2 / 8 



fc2 fc2^2 



{a + I3)r]~ ^{r + sf k^ ) w^fc^ 



— \/2l3r]ijj'^k'^ — —j=^ (r + s) + 2^*^^^ "^^^ 

- {sa - r/3) cj + (a - ufc^) (a + /3) ) u?k 
\2v2 V2 / 

(a + 2/3) ?7a;2fc2 (^s^u? + 2 (a - uk^f - A^ k^ 

-ii ((r + s)a + 2r/3) fc^w^ + (a - itfc^) /3k'^u!^ - ^iAk'^ (r + s) w 
i {{a + 3/3) auj^ + A (a + /3) fc^) a;^ 



2 7„2 



(A2a) 

(A2b) 

(A2c) 

(A2d) 

(A2e) 

(A2f) 
(A2g) 



Appendix B: Hofava-Lifshitz Gravity Projection Operators and Tensorial Relations 

In order to facilitate further use of the set of the degree-of-freedom and the parity preserving HLG projection 
operators we write out some useful identities satisfied by them. To avoid cluttering the notation, some indices of 
operators may be omitted. 



1. Vector field operators: N"^ — N"" 



Identities Among the Operators 



Tensorial Identities 



^ \^)ab- ^33 lab^-^ii V lab 



Sab =/'3r(ou+p^^(iu 

kakb = k^Pz3 (0)a5 



(Bl) 

(B2a) 
(B2b) 
(B2c) 
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2. Rank-2 Symmetric Field Operators, h — h 



Identities Among the Operators 

P''' (l)a.;c. = Pii (1++) + (1"") (B3a) 

P"" (2)a6;c. = (2++) + P2'2 (S^) (B3b) 

Tensorial Identities 

Sat,cd = I {SacSbd + SadStc) = P'"' (2) + P^"' (1) + P^"" (0^) + P^i" {0^) (B4a) 

SabScd = 2P/\'* (0'') + V2Pl'2^ {0"^) + V2P^I' (0"'') + P^2^ (0") (B4b) 

fcaA:fc<Scd + k.kdSab = (Fi'^ (0'") + ^2? (0"^)) + ^k^P^^^"^ (0") (B4c) 

fcaA^e^bd + kakdSbc + A:fcfc,<5„d + fc6fcrf<5„, = 2^ P'''' (1) + 4A;2p4'' (0") (B4d) 

kakbkckd = k^P^2^ (0") (B4e) 

{eecaSbd + eecbSad + SedaSbc + SedbSac) = {2P'^^ (3- + ) _ 2Plf (2+-) - P/'a'' + + P2l (l" + ))(B4f) 

{eecakbkd + Secbkakd + Eedafefefcc + £edbfcafcc) k" ^ 2k^ Vk^ {- Pif (l + -) + P^^^ (l- + )) (B4g) 

3. Graviton- Vector Operator, h — N 

Tensorial Identities 

P'"" i^ab:, - Plf (1^^) + (1-) (B5a) 

(5„fc/cc - %/2fc2pjWV (QS) ^ \/pP2'^f (0++) (B5b) 

I [kaS^b + kbS^a) = ^ (P^3^ (1++) + Pif {!-)) + ^P^i" (0++) (B5c) 

4. Graviton- Vector Operator, N — h 

Tensorial Identities 

P"" (l^).;.c = PsT (1++) + (1-) (B6a) 

kaSbc = \/2fc2P3^'' (0") + \/fc2P3^'' (O'^) (B6b) 

^ (kbSac + k,6ab) - ^ {P3T (1++) + (1--)) + VFPj^'^ (0-) (B6c) 
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